Runge-Kutta integration
https://www.youtube.com/watch?v=5CXhHx56COo